Z <- c(1, 1, 1, 0, 0,
0, 0, 0, 1, 1)
Y <- c(0.859, 1.930, 0.875, 2.944, -1.015,
-0.064, 1.624, -0.411, 1.048, -0.282)
dim_observed <- mean(Y[Z == 1]) - mean(Y[Z == 0])
se <-
sqrt(
var(Y[Z == 1])/length(Y[Z == 1]) +
var(Y[Z == 0])/length(Y[Z == 0])
)
2 * pnorm(abs(dim_observed), sd = se, lower.tail = FALSE)[1] 0.7382426